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We present a non-perturbative study of the (l+l)-dimensional massless Thirring model by using path integral 
methods. The regularization ambiguities -coming from the computation of the fermionic determinant- allow to 
find new solution types for the model. At quantum level the Ward identity for the IPI 2-point function for the 
' fermionic current separates such solutions in two phases or sectors, the first one has a local gauge symmetry that 

, is implemented at quantum level and the other one without this symmetry. The symmetric phase is a new solution 

^SJ ' which is unrelated to the previous studies of the model and, in the non-symmetric phase there are solutions that 

T-H ' for some values of the ambiguity parameter are related to well-known solutions of the model. We construct the 

, Schwinger-Dyson equations and the Ward identities. We make a detailed analysis of their UV divergence structure 

\l ' and, after, we perform a non-perturbative regularization and renormalization of the model. 
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^ 1 Introduction 

, The massless Thirring model is a well-known exactly solvable quantum field theory in (1+1) dimensions, as it has been 
1-^, shown, for the first time, by Thirring in his seminal paper pp. Thirring constructed the eigenstates of the Hamiltonian 
I and calculated some observable quantities and found them to be finite after renormalization. From that date, many 
studies and extensive investigation of the model have been done. Hagen (2) introduced an external field and found 
(— I \ the general solution of the Thirring model. Of course, the exact solution of the massless Thirring model means the 
possibility of an explicit evaluation of any fermionic correlation function, and it was carried out, for the first time, 
by Hagen ,2j, and afterwards by Klaiber 3 . In the Klaiber's analysis of the Thirring model, he found the operator 
solutions which are expressed in terms of a free massless Dirac field. On the other hand, Nakanishi ^ expressed the 
solution in terms of the free massless bosonic field and asserted that all Heisenberg operators should be expressed in 
terms of asymptotic fields from the standpoint of the general principles of quantum field theory. 

Recently a study has been made on the massless Thirring model has been made whose purpose was to investigate the 
equivalence between the massive Thirring model and the sine-Gordon model [Fj. On the other hand an extensive and 
detailed analysis of the model in its chiral symmetric phase was performed in [S] and its operator bosonic representation 
in 0. 

In previous works, the ambiguities have been observed in the fermionic current definition, as it was noted by 
Johnson 8 , Hagen 2^ and Klaiber , 3 , and it affects the coupling constant, therefore, in the massless Thirring model, 
the coupling constant is determined only when we define the current regularization. It was proved that the coupling 
constants are related to each other depending on the regularization. For example the coupling constant of the Schwinger 
definition j9j is given by 



X 
J3 
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where gj is the coupling constant of the Johnson definition ^ . Another study of ambiguity of the current regularization 
was performed in jl()| . These current ambiguities also appear in the massive Thirring model which still remains to be 
understood "11|. 

The aim of the present work is to show renormalized new solutions for the model obtained by using non-perturbative 
techniques via the functional integral method. The key to get such solutions are the regularization ambiguities that 
appear in the calculation of the fermionic determinant. Once the fermionic determinant is well-defined we are capable 
to solve the model without having to define the fermionic current [121 flSj . 
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Then, the regularization ambiguities aUow us to show that, at quantum level, there are two ways to get the full 
quantization of the masslcss Thirring model. The first one has solutions that are related to the various studies of 
the model just performed over the last years in which they do not consider the possibility of implementing a local 
gauge symmetry at the quantum level by exploring the ambiguities that appear in the calculation of the fermionic 
determinant or in the definition of the fermionic current. And the second one is a solution, as it will be shown, which 
results when we impose a local gauge symmetry at the quantum level. This imposition leads to an equation which 
fixes the ambiguity parameter as a simple function of the coupling constant. The unexpected result comes from the 
requirement of a local gauge symmetry which is not present at the classical level. 

The paper is displayed in the following way. In section 2, we present the model and shown the consequences of 
the regularization ambiguities in the solutions for the model. The section 3 and 4 shows a detailed analysis of the full 
quantization of the masslcss Thirring model in both phases mentioned in the previous paragraph. We do an extensive 
study about the calculus of the relevant Green's functions, the construction of the Schwinger-Dyson equations and 
Ward identities. Moreover, we make a careful analysis of the divergence structures and their regularization procedures, 
and finally we establish the renormalization and discuss their consequences. In section 5, we give our conclusion and 
perspectives. 

We use the following conventions: 

3™ = l = -.9", ."^ = l = -6oi, 
{7^7^} = 2g'^^ {7^,75} = 0, 75= 7V, P± = ^(l±75). 

2 The massless Thirring model 

The massless Thirring model ,1, is a theory of a self-coupled Dirac field defined by the following Lagrangian density 

= v;»w-| (v57^V')' , (1) 

where 5 7^ 0. The Lagrangian is invariant under the global chiral group U^,{^) x U^{1) at classical level, 

^ — >e'"V, — >e"^''^ip. (2) 
We quantize this fermionic model following the functional integral method, thus we write its generating functional 

as 

Z[ri,r]] = jvipVi) exp Jdx - | {i^l^i^f + # + i^-q^ . (3) 

The first step in the quantization procedure is to linearize the quadratic term in the fermionic current introducing an 
auxiliary vector field A^, such that the generating functional is expressed as 

Ziv, V] ^ J'DAf.Vi^Vip exp (^i Jdx + 4)tlj + ^A'^^ + f]ip + i'V^ , (4) 

where the presence of the massive term for the vector field precludes the existence (at classical level) of the following 
U{1) local gauge symmetry, 

^ _^ e'"V , V^e-'" , Af, ^ + d^a . (5) 

The fact that at classical level the model does not have a local gauge symmetry allows us to compute the fermionic 
determinant using a generalized regularization prescription ^3 1131 1141 1151 [T?)j to control the UV divergences. By 
following the point-splitting regularization implemented in |15j , we compute the fermionic determinant to be 



det(i^ + j^f — exp j dx 



a+1 d^^d" 
2 □ 



Au (6) 



where the parameter a characterizes the different regularization procedures used to control the UV divergences. When 
a = 1 the fermionic determinant is invariant under the local gauge transformation A^^ —t + 9^q;. 

In the following, we study the model for all the possible values of the ambiguity parameter a 7^ 1. 

After the fermionic integration the generating functional reads as 

Z[77, f]\= j VA^ exp (^iJdx^Af, (^^ + g^'' - A,^ exp (^-i J dx dy fi{x)G{x, y; A)r,{y)^ , (7) 



where G{x, y; A) is the Green's function of the Dirac's equation: (i^ + jjL) G{x, y; A) — 6{x — y), and it can be exactly 
computed, 

G{x,y;A)=expi-i dz Af,{z)j1{z, x,y)] P+Gpix ~ y) + exp i -i dz A^{z)j'^{z,x,y)] P-Gf{x - y) , (8) 



with Gf{x — y) being the Green's function of the massless Dirac's equation: ii^Gp^x ~ y) — 6(x — y); and the contact 
current jt^ (z, x, y) is 

j^{z,x,y) = {d^Td^)[DF{z-x)-DF(z-y)] , (9) 
where Dp(x) is the Green's function of the massless Klein-Gordon equation: DDf{x — y) = 5{x — y) . 
At quantum level the Ward identity for the IPI 2-point function of the fermionic current is given by 

q,T^''{q)^-{g+^\q'' . (10) 



a- 1^ 

then, we will have a local U (1) gauge symmetry if and only if 

a = l-^ (11) 
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such symmetry is reflected in the effective action Q which is invariant under the local U{1) gauge transformation 

A^^A^ + d^a . (12) 

The equation is important since the ambiguity becomes a simple function of the fundamental coupling of the 
model. The ambiguity that appears in the regularization process of the fermionic determinant or in the definition of 
the fermionic current is fixed by the requirement of a new local gauge symmetry at quantum level. Therefore, the 
model at quantum level displays two phases: one of them does not have the local U{1) gauge symmetry and the other 
one has. In the following we analyze both phases separately. 

3 The non— gauge invariant case: a 7^ 1 — y 

We begin the quantization process from the generating functional defined in and we include the external source 
for the fermionic current, thus, we have 

Z[77, fj, Gf,] = JvAf^V^Vi, exp (^i jdx i>{i^ + 4 + + ^^m^'' + # + ■ (13) 



In this situation we do not need a gauge fixing term for the vector field because its effective action Q is not gauge 
invariant under H12|) . Thus, considering the fermionic determinant in © for a 7^ 1 — — we can compute all the 
correlation functions and we can also construct the Schwinger-Dyson equations (SDE) and the Ward identities. 



3.1 The Green's functions 

3.1.1 The fermion propagator G{x — y) = {Q\T'4>{x)^l){y)\{)) 
The fermion propagator reads as 



ITT f dk 1 — e~^^'^^~y^ ' 
G{x -y)^^exp[ j ^ ]GAx-y), (14) 



where 



b=- + -^ (15) 

The fermion propagator has ultraviolet divergences and it will be shown that after its regularization and renor- 
malization procedures guarantee the full finiteness of the model. 



3.1.2 The current propagator J^^{x - y) ^ (0|Tj^(a;)jj,(?;)|0) 

In momentum space the propagator of the fermionic current = ^p^'^ip is calculated to be 

+ (16) 

9 + — T 9 + 7 

a + I a — I 

and it is free of UV divergences and, the tensors T^i, and i^j, are the transversal and longitudinal projectors defined 
as 

T,.^g,.~^ , and L,,^^. (17) 



3.1.3 The vertex function G^'{x,y]z) = {Q\Tip{x)il){y)j''{z)\Q) 
It is calculated to be 



In momentum space it reads as 



G^(p,-p-g;g)=z -— -— \ + -\75 Gip + q)-G{p) . (19) 

V 9ib - 1) r 9b J ^ -I 

where G{p) is the Fourier transforms of the fermionic propagator (|14|l . Thus, we can note that the UV divergences of 
this Green's function come from the fermionic propagator. 

It can be shown that 2n-point Green's functions for the pure fermionic current are UV finite, and all the 2n+m-point 
Green's functions involved 2n fermions and ni current insertions have logarithmic UV divergences. 

3.1.4 The IPI functions 

We write the IPI in momentum space. Thus, the IPI 2-point fermion function r(p) which is the inverse of the fermionic 
propagator is defined as 

V{p)^iG-\p). (20) 
The IPI fmiction 2-point function of the fermion current T^'^{q) is given by 

T^M) = I J-^{q) = - (.9 + ^) - [9 + ^) (21) 



r^{p,~p~q;q)^{ 1^+ l^j,] r{p + q)-T{p) (22) 



The IPI 3-point function r'^ (p, k; q) related to the vertex function G^{p, k; q) is 

' 27r 2tt q' 

a — 1 a + 1 q 

3.2 Schwinger— Dyson equations 

We start from the generating functional that after integration of the vector field can reads as 

Making the functional variation of the generating functional with respect to the fermionic field V' we get the 
quantum equation of motion 

where we have used the functional representations for the fields 

^ = "4 ' ^^^^ ' = ^''^ 



We perform a functional derivative with respect to r](jj) in the functional equation H24|) and setting the external sources 
to zero we obtain the SDE for the fermion propagator 



t^:,G{x - y) ^ iSix ^ y) + grG^{x,y;x) (26) 

which is satisfied by the Green's functions computed in 1)14(1 and H18|) from the generating functional when we 
consider the fermionic determinant 10) for a ^ 1 — 

As we mentioned above, the fermion propagator has UV divergences that are shown explicitly when we write the 
Schwinger-Dyson equation in momentum space 



Gip)--^-^J^fik)-JGip~k). (27) 



where the function f{k) reads for 



Iterating the SDE we generate an expansion in h or equivalcntly in powers of f{k) 



3 ^liSr^^'^^'^ri 'IwrWr^^'^^^'^l'^J^^^J^ 



which is a semi-perturbative expansion just as it was shown in |17j for the anomalous and chiral Schwinger model 
cases. A power counting analysis of the series clearly shows the existence of a logarithmic UV divergence in each order 
or loop, such property of the series guarantees that the model is a fully renormalizable field theory. 

3.3 The Ward identities 

We start from the generating functional 1(13(1 and we make the following fermionic local phase transformation 

^^e*^(")V , ^^^e-'^'^^K (30) 
In our framework the fermionic measure is not invariant under such transformation changing as 

V^pVt/j J[A, A] VtpV^P (31) 

where J[A, A] is the Jacobian given by 

det {id + 4 + (p) f.fj^^^^^A , ~ 1 \\2 



= det(.aV4+V-^A) = e^v(^Jdx^iA, + C,)^^X^^i^,Xr), (32) 
then, the generating functional in 1(13(1 reads as 

Z[rj,fj,Cf,] = JvA^Vi;Vij J[A] exp (^i Jdx i; {id + 4 + (p- ^X) i' + Af.A'' + fje'^^'=^ij + i;e-'^^'=^r]^ . (33) 
Next, we do the following transformation in the vector field — > + S^A and after integrating it we obtain 

Z[r], fj, C^] = jv^jVi) exp [i jdx ^i^^) - | (VW')^ + f^e'^'^^'H + ^e-'^^''\ + C^,j>^ (34) 



exp [ i dx 



a-1 6-1 , 



2tt 



6-1 (6-1) 
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By considering A(a;) as an infinitesimal function and expanding to first order in A, we get the functional equation 
satisfied by the generating functional 113() 

where we have used the functional representation ((25(1 for the fields. 

The generating functional of the connected Green's functions W[i],fj, C^] is defined as being W = —i\nZ and the 
generating functional of the IPI functions T[ip, , both are related via the following Legendre transformation 

T[,p,i;,j'^] = W[^,f},C^]- fdxiv^ + i^v + C^j'^), (36) 



and using the functional relations 

, 5W J 5W 5W _ 6T 5T ^ 5T 

^ = W ' ^ = ' '^^^ ' ""^H ' ' ^^^^ 

we obtain the functional equation satisfied by the generating functional r[?/;, j^] from which we get all Ward identities 



The first Ward identity is for the IPI 2-point function r'^'^{q) of the fermionic current which reads as 

q,T^-'iq) = -(^g + ^y- . (39) 

as it can be shown explicitly from the equation (|21|l . We can note that the current propagator is transverse for the 

critical value a ~ 1 and, as we mentioned in Section 2, it implies in the existence of a local U{1) gauge symmetry 

for this value of the ambiguity parameter. 

The second Ward identity which relates the IPI 3-point function r^{p,—p — q;q) and the IPI 2-point fermion 
function T{p) is given by 

q^T^^ip; q) - \t{p + q) ~ V{p)\ , (40) 

it can also be shown easily from the equation (|22() . This identity guarantees that all the divergences come from 
the fermion propagator. Thus, we can see that the ambiguity parameter does not renormalize and consequently the 
coupling constant either. 

3.4 Regularization and renormalization of the UV divergences 

We will study the UV divergences structure of the pure fermionic 2n-point Green's functions (correlation functions 
without current insertions) to identify their origin and later to perform their regularization. The (2n + m)-point 
Green's functions, where m is the number of current insertions, have similar divergence structure. Thus, as we will 
see, a simple renormalization of fermionic wave function is sufficient for all Green's functions becomes finite. 

Then, in the generating functional lO we decompose the vector field Af^ in their longitudinal and transversal 
components 

= d^P - M > (41) 
getting the following expression for the generating functional 

Z[n,f]] = jvpV(j) exp^i jdx ^^-^d^,pd^p- d f,(j)d^ (f^ exp (^i J dx dy 7]{x)G{x,y] p,(t))-q{y)^ , (42) 

where G(x, y; p, (j)) is given by 

G{x,y;p,(l))^^exp(^ijdz \p(z) ~ e(l){z)] j{z,x,y)^ P^Gpix ~ y) , (43) 

and the contact current j(z, x, y) = S{z — x) ~ 6{z — y). From the generating functional H42|l we can compute all pure 
fermionic Green's functions, for example, the fermionic propagator 

Gix-y) = if Vp exp U fdz ^^-^ df,pd''p + p{z)j{z,x,yyj (44) 



^ y !?(?!) exp (^-i jdz ^ 9^ 



09^0 + ec^{z)3{z, X, y) Gf{x - y) 



We can note that the functional integration on the fields p and (p are factorized. After performing the functional 
integration we can see that the ultraviolet divergences appear due to the bad behavior at high energy of the propagator 
of the transversal and longitudinal components of the vector field. The same situation happened when we compute 
all others higher order 2n-point and (2n + rn)-point Green's functions. 



As we have already known the origin of the ultraviolet divergences, our next step will be to construct a non- 
perturbative scheme of regularization inside the generating functional. Thus, we begin the regularization process from 
the generating functional (|13|) setting = which reads 



Z[77, fj] = JvAf^VipVifj exp (^i Jdx {id + 4) i} + ^ A^,A^' + f]ip + i/"7^ • (45) 

We decompose the vector field as in (|41ll and make the following fcrmionic transformation 

^^eV-*750^ ^ ^ ^ ^ f^-w-ii^4> _ (46) 

As it is well-known, the fermionic measure is not invariant under a chiral transformation |18) . On the other hand, 
we can also show that it is not invariant under a phase transformation |14l IT!?] due to the generalized prescription 
used to compute the fermionic determinant. Thus, the fermionic measure under the transformation (|46|) changes as 
'D'ip'D'tp J[p, (j)]'DipT>'ilj , where J[p, cj)] is the Jacobian given by 

Inserting 141() . H46() and (|47|l into the generating functional (|45() . we get 

Z[r], fj] = jvpV4)V%l;Vi} exp jdx ^i^^j + df.pd'^p ~ df^cpOf^cj) + fie'P~'''"^^tp + i/ie"*''-*'^^'^?]^ . (48) 

The fundamental observation here, is that we could make everything finite if we had a better UV behavior for the p 
and (j> propagators. Then, by following the scheme showed in [T?) . we add to the generating functional two regulator 
fields a and f3 with large masses A^ and A^^ respectively 120] (A^ +oo). Also, the interacting terms of p and (p fields 
are changed as p ^ p + a and (/> — > + /5. In this way, we get to define the regularized generating functional as being 

z'^[T],f]] = JvpVcpVijViljVaVp exp(^i Jdx-^^^,<j)^^'(j)+^ df,f3d^'(3-Al(3^^ (49) 

exp (^t Jdx ^ d^pd'^p - ^ [^^a^^'a - Ala^]^ 

exp (^i jdx^pi^Tp + j7e'[p+"~T5(0+/3)]^ _^ ^^-^[p+o^+l^M+P)]^ _ 

Now, before carrying out the integration of a and (3 fields, we make a change of variables p ^ p — a and cj) ^ cj) — p. 
Thus, we obtain the following expression for the regularized generating functional 

(50) 



e^P /'^^ ^ + ^> - ^ + ■ 

In this way we were able to improve the ultraviolet behavior of both propagators. Being able to regularize those 
propagators, we are going to express the regularized generating functional in function of the original fields as in l|45|) . 
Doing the change of fermionic variables inverse to that one made in (|46|l . and expressing fields p and (f> as functions of 
the field by using H41|l . we obtain 

Z*[r?,7y]= JVA^V^PV^ eicp (^i Jdx {^{1^ + 4)4> + ^A^A^ + ^ {O- Af - {d- Af + + ^tj^ (51) 

Comparing the Lagrangian densities that appear in the generating functionals H45|l and H51|l . we can see that the 
regularization process has given rising to two new terms which depend explicitly on the Pauli-Villars cut-offs, such 
terms are enough to turn completely finite the model. We can to compute all the regularized fermionic Green's 
functions starting from H50|l or (|51() . 

By integrating the auxiliary field A^, we get the regularized generating functional in terms of the fermion fields 
alone, 

[V, V] = Vi^Vi) exp(i dx Vii^V - | (V'T^V')^ + | i^i'^i'Dl^iiYi^ + # + ^/;?7 J , (52) 



where the operator D carries all the information about the non-local counter-terms which regularize the fermionic 
theory, 



ttA^ 



9{b-l) 



7rA2 ' 



(53) 



as we can see, it shows explicitly the longitudinal and transversal contributions that regularize the respective compo- 
nents of the fermionic current -current interaction. 

3.4.1 Regularized and renormalized fermion propagator 



We start from the regularized generating functional 1)51) |l. integrating over the fermionic fields it reads as 



r]{x)G{x,y;p,4>)r]{y) , 



Z^hv] = JvpVcj) exp (^i jdx </.□(□ + A2)0 - ^-i- pU{U -f A^)^^ exp [-i Jdxdy 

(54) 

where the Green's function G{x,y; p, (p) is given by (|43() . Then, the regularized fermion propagator in configuration 
space is given by 



G (x — y) — i exp [ i 



dk 
(2^ 



7rA2 



1 



1 



6-1 fc2(fc2-A2) 6 F(fc2_A2) 



Gf{x - y) , 



(55) 



performing the momentum integration in the exponential and introducing the arbitrary mass scale M , the propagator 
is 

1 1 

1 



G {x,y)^ 



1 



2TTi {x-yf-i£ 
where the parameter c is defined as being 



4(6-1) / A|\46 
i M2 ] 



- A'Pix-yf + i£ 



(56) 



c = 46(6- 1) = ( y + a ) - 1, 



(57) 



we can see that c G [—1, -f oo > and c ^ € < — oo, —1] U < 0, -l-oo >. 

In this way, we can renormalize exactly in configuration space [21|.The renormalization procedure establishes the 
following relation between renormalized and bare (regularized) IPI Green's function, T = Z^, F , where Z^p is the 
fermionic wave function renormalization constant that relates the bare and renormalized fermionic fields, V'fcare = 
Z^fjipji. Thus, from (|56|) . we choose the fermionic wave function renormalization constant to be 



AM 4(6-1) / An 46 



M2 



(58) 



It can be shown that a renormalization of the fermionic wave function is sufficient to get renormalize all the Green's 
functions of the model. 

Then, the renormalized propagator of massless Thirring model in configuration space is read as 



G (x-y)^ 



i2TT {x — 2/)2 — ie 



-l/c 



(59) 



in the limit g we also obtain the propagator of a free massless fermion field. 

Then, by computing the Fourier transform we get the renormalized fermion propagator in momentum space, 

1 



n I (AP\ c .,;^/, F(1-1/c) 



(60) 



where M — 2M. The dynamical or anomalous dimension of the massless Thirring fermion field is l/c. By choosing 
an adequate value for the ambiguity parameter we can reproduce results such as has been obtained in [Tl 1221 1^ I24| . 



This value is got when we set c = '2.'^ / 9^ thus, the ambiguity parameter which reproduces the result to the anomalous 
dimension found in the literature is 

27r I 27r , , 

..:=--±^l + -. (61) 

In this value we will obtain a result which agrees well with Glaser's analysis of the massless Thirring model |22II23| 
and spectral analysis of two-point correlation functions in 1+1-dimensional quantum field theories developed by 
Schroer 

The renormalized propagator 161) II has singularities in - = ±1, ±2, . . .. And, if we set c = 27r/g the propagator 
(|60|l has singularities in g — ±27r such as it was shown by Klaiber analysis p . 

4 The gauge invariant case: a = 1 — y 

As we have observed in Section 2, the Ward identity for the IPI 2-point function of the fermionic current reveals the 
existence of a local U{1) gauge symmetry for a = 1 — For this value such symmetry is explicit in the effective action 
for the vector field in the generating functional iQ, i.e., we can not to perform the functional integration. Therefore, 
we need to construct a well-defined generating functional for this gauge theory. Such construction is making via the 

Faddeev-Popov procedure |25| with Lorentz gauge fixing condition f{A^) — —= d- A in the generating functional 
(0J. Thus, the resulting generating functional for this gauge invariant theory reads as 

Z' [?7, 77, C^] = j VA^Vi^Vi; exp (^i J dx ^{i^ + 4 + (^yj + A^.A^' - ^ {d- A)^ + f]^p + ^tj^ , (62) 

where we have introduce the external source for the fermionic current. From this generating functional we can 
compute all the Green's functions and we can also to construct the Schwinger-Dyson equations (SDE) and the Ward 
identities. 

4.1 The Green's functions 

4.1.1 The fermion propagator G"" {x — y) = {0\T'ip{x)'ip{y)\0) 

We set = in the generating functional H62|l and integrating the fermionic degrees of freedom we obtain 

Z'[r^,fj] = jvA^ exp jdx \a^ ^ g'^" - ^ ^ + ^ d'^d'' A^^ exp (^-i Jdxdy f]{x)G{x,y; A)r^{y)^ , (63) 

where the Green's function G{x,y;A) is given in Q. From this, we calculate the fermion propagator which becomes 

G^(x-y)=* exp 1^(1 + [l - e^^^-^)] ) G^x - y) , (64) 

We can see the existence of logarithmic UV divergences coming from the term in the integral proportional to k~'^ 
and, we also observe that the term proportional to the gauge parameter ^ generates serious infrared divergences which 
can not be renormalized. These infrared divergences proportional to ^ are also observed in the all Green's functions 
computed from the generating functional 162|l . then, we use a gauge free of infrared divergences ^ — 0. 

4.1.2 The current propagator J^^ — {0\Tj^{x)j,^{y)\0) 

In momentum space the propagator of the fermionic current is calculated to be 

-(i--)t. + -(^ 

9 \ 9 ) ^ 9 \9k 

it is free of UV divergences. 



iJ^Al) = 7 ( 1 - T ) + 7 ( - 1 ) ^M.- (65) 



4.1.3 The vertex function G^(x,y;z) = {0\Ti}{x)'tp{y)jf,{z)\0) 
The vertex function is calculated to be 



in momentum space it reads as 



GJp, -p-q:q) =i 



gq J q g q 



G\p + q)-G\p) 



(67) 



where (p) is the Fourier transforms of the fermionic propagator H64|) . Thus, we can note that the UV divergences 
of this Green's function come from the fermionic propagator. 

4.1.4 The IPI functions 



In momentum space the IPI 2-point fermion function F {p) is 

T'ip)^i\G'{p) 



Also, the IPI 2-point function for the fermionic current T (q) is 



T 



(68) 
(69) 



- 1 



gq 



The IPI 3-point function, F^(p, k; q) related to the vertex function G'^Ap, k; q) is given by 



r^(p, -p-q;q)= \-9^ + ^ 



(70) 



4.2 The Schwinger-Dyson equations 

We start from the well-defined gauge fixed generating functional (|62|) which after the integration of the vector field 
reads as 



Z' [ry, fj, G^] = V^V4, exp (i dx ^{i^ + (P)^ - ^ it^f + f fOU'^ + f^i, + 4,7^ 



2" 

where the third term appears due to the gauge fixing condition, and the operator is given as 



□ + i 



(71) 



(72) 



By using the functional representation for the fields H25|l in the equation of motion for the fermionic field, we get 
the quantum equation of motion for the generating functional (|71|l 



6n{x) ''^^ 5f){x)5Gf^{x) 57){x) ' ^' 6f]{x)^ SG,{x) ^ 

The last term makes a difference between the SDE's (|73|l and (124(1 but it arises due to the introduction of a gauge 
fixing term into the generating functional (|62|) for well-defined it. Thus, in a formal treatment the gauge fixing term 
is not present and therefore the SDK for the gauge invariant phase is exactly the same as in (I24|l for the non-gauge 
invariant phase . 

Thus, the SDE for the fermion propagator in the gauge invariant phase is given by 

"^''r — ii\ = 1AI T — ii\ 4- n^f( , ,, „ „, „, ,^ , 



zC^(x)7'^- 



Z [,7,r7,C"^]=0 (73) 



i^xG\x-y) = i6{x - y) g-^''G^ {x,y;x) - g^^' Wva Gl^{x,y]z) 

Z—^X 



(74) 



Both the fermionic propagator (|64ll and the vertex function 166(1 computed from the generating functional ((62(1 are 
also solutions for the quantum equation of motion of the massless Thirring model. 
In momentum space the SDE for the fermion propagator reads as 



where the function / (k) reads for 



(75) 



(76) 



Similarly to the case non-gauge invariant ((29(1 by iterating the SDE it generates an expansion in h or equivalently 
in powers of (fc) . The power counting analysis of the series shows clearly the existence of the same logarithmic UV 
divergences in each order or loop what indicates that the model is a fully renormalizable field theory. 



4.3 The Ward identities 

We start from the well-defined generating functional Ht)2|) and we make the following fermionic transformation 



(77) 



one more time it is worthwhile to note that in our framework the fermionic measure is not invariant under this 
transformation changing as VipVip J [A, A] VipVip where J [A, A] is the Jacobian given by 



y [A, A] = exp ^ jdx - ^ [A 



(78) 



The generating functional (|62|l after the fermionic transformation reads as 

[rj, f], Cf,] = / VA^V^Vijf [A] exp ( i [dx {id + 4 + (p- ^X) + ^ AA - ^ {dA) 



2g 



2C 



At once we make the transformation A^ A^ + 9^A and after integrating over the vector field we obtain 



[77, fj, = IV^Vi, exp ( I jdx ^{i^ + (p)^-'- (V'7''^)' + ? fDl.f + jye'^^^V + 



iX{x 



(79) 

(80) 



exp 



/ 



v 



dx-- C^a^A + A— 
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1 



\ 



9 ) 



Considering A(a;) as an infinitesimal function and expanding the exponential to first order in A we get the functional 
equation satisfied by the generating functional 



/ 



V 9 



Z^ [rj,fj,C^]^0 



(81) 



where we have used the functional representations for the fields (|25|1 . 

We define the generating functional of the connected Green's functions [r/, 77, C^] as being w"" ~ —i\nZ^ and 
the generating functional of the IPI functions ip, j^^] which are related via a Legendre transformation similar 

to H36(l . Using the functional relations given in (|37|l . we obtain the functional equation satisfied by the generating 
functional [ip, ijj, j'^] from which we get all Ward identities 



dipix) 



dipix) 



9 ^dj^{x) 



9 



d^^rix)^o. 



(82) 



If we observe the same functional equation in the non-gauge invariant phase H38|l , the equation above can be obtained 
27r 

by setting a = 1 . The last term appears due to the gauge fixing condition used to define the generating functional 

9 

of the Green's functions in this phase. 

The first Ward identity is to the IPI 2-point function T^n,^{q) for the fermionic current which reads as 



q<M^9 — 



(83) 



as it can be shown explicitly from the equation H69|l . Thus, the right-handed side above is gauge dependent what 

27r . . __ 

guarantees that the current propagator is transverse for the critical value a = 1 , as it was noted in (|39|l . 

9 

The second Ward identity relates the IPI 3-point function r^{p-,q) and the IPI 2-point fermion function F (p) is 
given by 

q^^rl [p- q)^-g [f' {p + q) - (p)] (84) 

likewise it can also be shown easily from the equation (|70|l . This identity shows that the coupling constant does not 
renormalize due to a renormalization of the fermionic propagator guarantees the rcnormalization of this IPI 3-point 
function. 



4.4 Regularization and renormalization of the UV divergences 



We will study the divergences of the pure fermionic Green's functions (Green's functions without current insertions) 
following the procedure achieved in the subsection 13.41 for the non-gauge invariant phase. Therefore, starting from 
the generating functional H63|l and decomposing the vectorial field A^^ in their longitudinal and transversal parts (|41|l 
we arrive to the following expression for the generating functional for the all pure fermionic Green's functions 



Z'[7^,fj] = JvpVq^ exp (^ijdx - ^ d^ct^df^cf)^ ^ {Dpf^ exp (^-i jdxdy n{x)G{x,y- p,4>)il{y)^ 



(85) 



where the Green's function G(x, y; p, <j)) is given by equation (|43|l . Thus, in order to compute the fermionic propagator 
we take functional derivatives with respect to the external fermionic sources and setting r] = fj = we obtain 



G\x-y) = ijvpV(^ exp(^zjdx ^ 1- d^c/^d'^^ - (apf^ G{x,y;p,, 



(86) 



where the integrations over the field p and (p are factorized. It can easily be observed that the UV divergences come 
from the integration on the transversal field (j) due to the bad behavior of its propagator at high energy, thus, we need 
to improve its behavior at high energy to be able to control the UV divergences. 

Therefore, we start the regularization process from the generating functional 1)62(1 setting — 0, then, we decom- 
pose the vector field as in 141|) and make the fermionic transformation H4ti|) which changes the fermionic measure as 
VipV^ [p, (flVipVij; , where [p, cf] is the Jacobian given by 

f [p, 0] = exp (^ijdx - j- d^pd^p - ^ (f - l) a^-^a^^) . (87) 

Finally, we get the following expression for the generating functional 

Z\ri,f]] = JvpV<j)VipV^P exp^i Jdx V^i^V" ^ (Op)^ - df,(f>d^'<j) + f]e'''-'^'''''ip + i;e~'P-''^'"''r]^ 

Then, we following the regularization procedure described between the equations and -for the non-gauge 
invariant phase- but improving only the 0-propagator. Such procedure leads to the following regularized generating 
functional for the gauge invariant phase 

Z^^'lr^^fj] = JvpV(l)VijV^ exp ^i Jdx ^i^^j- ^ (Dp)^ + (f>D{D + A^),/) + T^e'P-'^'-'^'ip + i^e-'P-''^"'^!^ . (89) 

Such regularized generating functional expressed in terms of the original fields reads as 

I . ... 1 ..T 1 

or, we can also express it as a pure fermionic theory integrating the auxiliary field 



Z "' [r,, fj] = jVA^V^V^j exp { i j dx ^^{i^ + 4)V + ^^m^"" + 1^ (9- A'f - — (5- Af + fjib + ^r)\ (90) 



where the operator is defined as 



^f.u - T- + — r2— • (y^J 
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The third term in the regularized generating functional H91|l has two contributions: the longitudinal one which has 
already observed in (|7TT) is due to the gauge fixing condition and the transversal one that arises due to the regularization 
procedure of the UV divergences. Thus, the gauge fixing condition acts as a natural regulator of the longitudinal part 
of the current-current interaction and the Pauli-Villars method was necessary to regularize the transversal component 
of the interaction. 



4.4.1 Regularized and renormalized fermion propagator 



From the regularized generating functional H89|) we can compute all the regularized pure fermionic Green's functions 
of the local gauge invariant model. Thus, for arbitrary gauge ^ the regularized fermion propagator is given by 



g"\x-v) 



I exp 



dk 
(2^ 



1 

fc4 



Gf{x - y) , 



(93) 



Now, going to the ^ = gauge we compute the remaining integration to obtain the following expression for the 
regularized fermion propagator 



G {x-y}^ 



27ri {x — y^ 



IS 



Af2 



1 



(94) 



where we have introduced an arbitrary mass scale M . 

Immediately, we renormalize in configuration space by choosing the fermionic wave function renormalization con- 
stant as being 



— 



(95) 



Thus, in configuration space the renormalized propagator of massless Thirring model in, its gauge invariant phase, 
reads as 



G (x-y) 



( ~M'^{x-yf + i^ 



i27r (x — y)"^ — ie 

Computing the Fourier transform we obtain renormalized fermion propagator as being 



(96) 



1 



t/4 



£(5/4) 

r(3/4) 



(97) 



This propagator deserves some comments. It does not have any explicit dependence in the coupling constant nor in 
the ambiguity parameter. Moreover, it is similar to the infrared behavior of the fermionic propagator of the Schwinger 
model j2f)l 1271 1^. which has the confinement phenomenon as one of their characteristics. And, the Schwinger model 
presents chiral symmetry breaking due to the correlation functions of the chiral densities break the cluster property 
which leads to 0— vacuum structure. Thus, two important quantum field theoretical fenomena could also be natural 
characteristics of the massless Thirring model in its gauge invariant phase. 



5 Remarks and conclusions 

The massless Thirring model was analyzed in a non-perturbative way by using functional integral method. The key 
of the present analysis is the calculus of the fermionic determinant which was made using generalized prescription 
to control the UV divergences that appear in its computation. This arbitrary procedure introduces an ambiguity or 
Jackiw-Rajaraman parameter a and it allows to obtain the new solutions for the massless Thirring model. 

Next, we explain the contribution of the ambiguity for the new solutions: the fermionic determinant is defined by 
the following functional integral 

det {i^ + 4) = y Vil;V4> exp (^i Jdx 4> {i^ + 4) ij^j (98) 

if we make the following local phase transformation 

^^e*'''="V , ^^i^e-'P^''^ (99) 

in our prescription the fermionic measure is not invariant under such transformation and there is a nontrivial Jacobian 
given bv pilTl] 

det(i^ + 4) f ■ f , o, — ^ A„ a — 1 ^ ^„ A 

= det(^^ + 4^P) = VJ'' ^ ^''^''^ "-^''^j ■ ^'''^ 



It shows the exphcit breaks of a local U{1) symmetry which, in general, is associated to the conservation of the 
fermionic current. But such one local symmetry does not exist for the massless Thirring model at classical level, 
thus, we are free to compute the fermionic determinant in an arbitrary prescription characterized by the ambiguity 
parameter a. 

At quantum level the Ward identity for the IPI 2— point for the fermonic current, for a ^ 1, allows to choose a 

value for the ambiguity parameter in which it is possible to implement a local C/(l) gauge symmetry for the model. 

This fact is reflected after performing the fermionic integration, at this point it arises a problem when we wish to 

perform the functional integration over the auxiliary vector field for all values of the ambiguity parameter due to 

27r . . . 

there is a critical value a = 1 where it is not possible to be performed. For this value, the effective action 

9 

of the vector field gains a local U{1) gauge symmetry and we use the Faddeev-Popov technique to well-define the 

27r 

functional integration. Consequently, we named as non-gauge invariant phase to the solutions with a ^ 1 and 

9 

27r 

gauge invariant phase to the solution with a = 1 . We also show that the Green's fimctions. in both phases, 

are solutions of the Schwinger-Dyson equation for the massless Thirring model which guarantees the validity of our 
approach. 

We identify clearly the origin of the UV logarithmic divergences present in the Green functions and make a 
non-perturbative regularization into their generating functional via the Pauli-Villars method. Such detailed analysis 
of the divergence structure is fundamental to carry out the covariant non-perturbative regularization of the the 
UV divergences, as it was observed in The UV divergences has a non-perturbative character due to they can 

not be obtained via the usual perturbative technique which uses the coupling constant as parameter expansion; this 
divergence type has already observed in the anomalous and chiral Schwinger model |19| . Accordingly, the regularization 
procedure modifies the components transversal and longitudinal of the current-current interaction yielding us the non- 
perturbative counter-terms which make possible the full renormalization of the model in both phases. 

A straightforward calculation for the higher fermionic correlation functions from the regularized Lagrangian (|52|l 
and H91|l shows that they have already regularized. This fact is notably important for the full renormalization of 
the model, because it implies in the non existence of additional divergences in the theory. Thus, in both phases, the 
coupling constant g does not renormalize such as it can be shown from the Ward's identities (|40|l and H84|l what was 
also observed in the literature The Green's functions for the composite operators, as for example, chiral densities 
need an additional renormalization procedure what will be reported in a following paper. 

In the non-gauge invariant phase there are some values of the ambiguity parameter that reproduce the solutions 
found in the literature computed by different methods or techniques. On the other hand, the gauge invariant phase is 
a solution which is implemented, for the first time, in the present work and it can not be reproduced for no value of 
the ambiguity parameter or coupling constant related to the non symmetric phase. 

In the gauge invariant phase, we note that for a general value of the gauge parameter ^, the Green's functions are 
plagued of infrared divergences but in the ^ = gauge they are infrared finite, then, we can work in this gauge. On 
the other hand, the non-symmetric phase is free of infrared divergences. 

The gauge invariant phase shows some peculiar features, the first feature is a simple equation Hll() which relates 
the ambiguity parameter a and the coupling constant g. Possibly, it could be fixed the ambiguity problem in the 
definition of the fermionic current in the present model or given a guide to solve it in the massive Thirring model 
This possibility comes at quantum level where the model gets a local gauge symmetry which allows to fix 
the ambiguity related to regularization freedom originated from the classical level. The other one is the fermionic 
propagator structure which resembles the behavior of the infrared limit for the fermion propagator of the Schwinger 
model j26[ 1271 1^ . As it is well-known the Schwinger model displays many characteristics of (1+3) dimensional more 
realistic gauge theories, such as confinement or 0-vacuum structure, such characteristics could emerge in the present 
phase. The study of the Hilbert space structure of the symmetric phase would allow understanding if there is fermion 
confinement and would identify the physical operators establishing clearly which symmetries are present and which 
are violated at quantum level. These questions mentioned above are at the moment under investigation and results 
will be reported elsewhere. 
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